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We analyze a simple model of deterministic diffusion. The model consists of a one-dimensional 
periodic array of scatterers in which point particles move from cell to cell as defined by a piecewise 
linear map. The microscopic chaotic scattering process of the map can be changed by a control 
parameter. This induces a parameter dependence for the macroscopic diffusion coefficient. We 
calculate the diffusion coefficent and the largest eigenmodes of the system by using Markov partitions 
and by solving the eigenvalue problems of respective topological transition matrices. For different 
boundary conditions we find that the largest eigenmodes of the map match to the ones of the simple 
phenomenological diffusion equation. Our main result is that the difffusion coefficient exhibits a 
fractal structure by varying the system parameter. To understand the origin of this fractal structure, 
we give qualitative and quantitative arguments. These arguments relate the sequence of oscillations 
in the strength of the parameter-dependent diffusion coefficient to the microscopic coupling of the 
f~>) , single scatterers which changes by varying the control parameter. 
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PACS numbers: 02.50.-r,05.40.+j,05.45.+b,05.60.+w,47.52.+j,47.53.+n 



I. INTRODUCTION 



Over the past years there has been a rapidly growing interest in trying to understand the mechanism of nonequilib- 
rium transport on the basis of dynamical systems theory One line of work is related to computer simulations, 

where interacting many-particle systems under nonequilibrium conditions, like shear or an external field, are in the 
center of the investigations [p~0|— p^4|| . Another line of research focuses on low-dimensional models like the random [l5|-|l^1 
or periodic Jl9|-pT[| Lorentz gas. An even simpler model which shares certain properties of the periodic Lorentz gas 
are two-dimensional multibaker maps |32]-|39t]. Lorentz gases and multibaker maps have become standard models in 
the field of chaos and transport, since, on the one hand, they catch the physical essence of certain real nonequilibrium 
processes, but, on the other hand, they are still simple enough such that they can be analyzed in detail theoretically. 
The final point along this line of reduction of complexity is the problem of deterministic diffusion in one-dimensional 
q , maps, as it has first been studied by Grossmann and Fujisaka by Geisel et al. E3-Ej|, and by Schell, Fraser 



and Kapr al |4ql. Later, methods of cycle expansion [47H5jJ and other techniques have been applied to this and related 
questions [[52|-|56|| . The most interesting problem in this context is to describe the dependence of deterministic diffusion 
on varying a single parameter value. However, for computing parameter-dependent diffusion coefficients so far either 
the deterministic dynamics could not be treated in full detail, and thus the results were approximative, or the task 
. . 1 could be performed exactly only for a few simple cases of parameter values. 

The idea of this article is to apply methods of dynamical systems theory, as discussed by Gaspard et al. |2j||32|,[37j , 
to the problem of parameter-dependent deterministic diffusion in one-dimensional piecewise linear maps. In the 
remaining part of this introductory section, we will explain the term deterministic diffusion, and we will define the 
class of dynamical system we want to analyze. In Section |J, some important background of our approach will be 



discussed briefly, and in Section III , a method will be presented which enables the exact computation of deterministic 



diffusion coefficients for a broad range of parameter values. The result for the diffusion coefficient of the simple map 



considered here turns out to be surprisingly complex so that additional investigations, performed in Section IV, are 
required to understand the origin of this unexpected non-trivial diffusive behavior. 

This article is based on the work of Ref. plj], Chapter 2, a brief summary has been published in Ref. p9|. 



A. What is deterministic diffusion? 

In a first approach, one may think about diffusion as a simple non-correlated random walk J60| — 163[] , as sketched in 
Fig. [|: One starts by choosing an initial position xq of a point particle on the real line, e.g., xq = 0, as shown in the 
figure. The dynamics of the particle is then determined by fixing a probability density p(s) such that p(s) = p(s) ds 
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gives the transition probability p(s) for the particle to travel from its old position xo over a distance s = x\ — xq to its 
new position x\. The same procedure can be applied to any further change in the position of the particle. For sake 
of simplicity, it is assumed that p(s) is symmetric with respect to s = 0. The quantity of interest is here the orbit of 
the moving particle, which is represented by a collection of points {xq, x\, X2, ■ ■ • , x n } on the real line, n refers to the 
discrete time, which is given by the number of iterations, or "jumps" , of the particle. 

For performing a simple non-correlated random walk, it is assumed that the probability density p(s) is independent 
of the position x n of the particle and of the discrete time n, i.e., it is fixed in time and space. In other words, 
there is no history in the dynamics of the particle, that is, the single iteration steps are statistically independent 
from each other. This is called a Markov process in statistical physics Such a random walk can be used 

as a simple model for diffusion: One starts with an ensemble of particles at, e.g., xq = 0, and applies the iteration 
procedure introduced above to each particle. If one separates the real line into a number of subintervals of size Ax, 
the number n p of particles in any such small interval after n iterations can be counted. Divided by the total number 
Np of particles, these results determine the probability density p n (x) for a sufficiently large number of particles. This 
quantity can be interpreted as the probability to find one particle at a displacement x after n iterations. Thus, p n (x) 
provides information about the macroscopic distribution of an ensemble of particles, or of the average position of one 
particle, respectively. The diffusion coefficient D can now be obtained from the second moment of the probability 
density p n (x), 

D= lim ±f-± , (1) 

where < . . . >:= J dx p n {x) . . . represents the probability density average. If the probability densities are Gaussians, 
the diffusion coefficient of Eq. (|l|) is well-defined. This is the case for a stochastic diffusion process as the one modeled 
here by a simple random walk. 

In contrast to this traditional picture of diffusion as an uncorrelated random walk, the theory of dynamical systems 
makes it possible to treat diffusion as a deterministic dynamical process: Here, the orbit of a point particle with initial 
condition xq may be generated by a chaotic dynamical system 

x n+ i = M(x n ) . (2) 

M(x) is a one-dimensional map which determines how a particle gets mapped from position x n to position x n +i, as 
will be introduced in detail below. The map M(x) plays now the role of the model-inherent probability density p{s) of 
the previous random walk: Defining M{x) together with Eq. (Q) gives the full microscopic equations of motion of the 
system. Thus, the decisive new fact which distinguishes this dynamical process from the one of a simple uncorrelated 
random walk is that here is uniquely determined by x n rather than having a distribution of x n +i for a given 

x n . This way, the complete history of the particle is taken into account. If the resulting macroscopic process of an 
ensemble of particles, governed by a deterministic dynamical system like map M (x), turns out to obey a law like Eq. 
(Q), i.e., if a diffusion coefficient exists for the system, this process is denoted as deterministic diffusion ]40|-|5l"[ . 



B. The deterministic model 

Fig. H shows the model which shall be studied in the following. This model has apparently first been introduced 
by Grossmann and Fujisaka fdlijjilj] . It depicts a "chain of boxes" which continues periodically in both directions to 
infinity, and the orbit of a moving point particle. Let 

M a : R -> K , x n 1 ► M a (x n ) = x n+1 , a > I , x n £ K , n E N (3) 

be a map modelling the chain of boxes introduced above, i.e., a periodic continuation of discrete one- dimensional 
piecewise linear fisty expanding maps with uniform slope. The index a denotes a control parameter, which is the 
absolute value of the slope of the map, x n is the position of a point particle, and n labels the discrete time. Since 
the map is expanding, i.e. a > I, its Lyapunov exponent In a is greater than zero. Thus, M a (x) is dynamically 
instable and may in this sense be called chaotic |36|. In order for the map to be chaotic and piecewise linear, it 
cannot be monotonic, so there must be points of discontinuity and/or non-differentiability. The term "chain" in the 
characterization of M a (x) can be made mathematically more precise as a lift of degree one, 



M a (x + I) = M a (x) + 1 , (4) 



for which the acronym old has been introduced |67[-|69|. This means that M a (x) is to a certain extent translational 
invariant. Being old, the full map M a (x) is generated by the map of one box, e.g., on the unit interval < x < 1, 
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which will be referred to as the box map. It shall be assumed that the graph of this box map is point symmetric 
with respect to the center of the box at (x,y) — (0.5,0.5). This implies that the graph of the full map M a (x) is 
anti-symmetric with respect to x = 0, 



M a (x) = -M a {-x) 



(5) 



so that there is no "drift" in the chain of boxes. For convenience, the class of maps defined by Eqs. (|), (|), and (§) 
shall be denoted as class V (where V stands for piecewise linear), and maps which fulfill the requirements of class V 
shall be referred to as class V -maps. In Fig. |2|, which contains a section of a simple class P-map, the box map has 
been chosen to 



M a (x) 



ax 



ax 

— a 



< x < \ 
\<x<\ 



a>2 



(6) 



cf. Refs. g57|j6ej]. This example can be best classified as a Lorenz map with escape [0-0. The chaotic dynamics 
of these maps is generated by a "stretch-split- merge" -mechanism for a density of points on the real line p% . As a 
class P-map, Eq.(p|), together with Eqs.([|), (Q), and (||), will be referred to as map C. Other class P-maps have been 



77,7; 



considered in Refs. p0| , ^l|, 47 57 7^j7^], another example is discussed in Refs. 

It has been proposed [^ | to look at the dynamics in this chain of boxes in analogy to the process of Brownian motion 
3211 : If a particle stays in a box for a few iterations, its internal box motion is supposed to get randomized and may 



resemble the microscopic fluctuations of a Brownian particle, whereas its external jumps between the boxes could be 
interpreted as sudden "kicks" the particle suffers by some strong collision. This suggests that "jumps between boxes" 
contribute most to the actual value of the diffusion coefficient. Brownian motion is usually described in statistical 
physics by introducing some stochasticity into the equations which model a diffusion process. The main advantage 
of the simple model discussed here is that diffusion can be treated by taking the full dynamics of the system into 
account, i.e., the complete orbit of the moving particle is considered, without any additional approximations. This is 
another way to understand the notion of deterministic diffusion in contrast to diffusion as obtained from stochastic 
approaches. That is, in the purely deterministic case the orbit of the particle is immediately fixed by determining its 
initial condition. 

One should note that the strength of diffusion, and therefore the magnitude of the diffusion coefficient, are related to 
the probability of the particle to escape out of a box, i.e., to perform a jump into another box. This escape probability, 
however, as well as the average distance a particle travels by performing such a jump, changes by varying the system 
parameter. The problem which will be solved in the following is to develop a general method for computing parameter- 
dependent diffusion coefficients D(a) for class P-maps. Here, map C will serve as a simple example. However, the 
methods to be presented should work as well for any other class V-map, supposedly with analogous results (see, e.g., 



II. FIRST PASSAGE METHOD 

The methodology of first passage, as it has been developed in the framework of statistical physics deals 
with the calculation of decay- or escape rates for ensembles of statistical systems with certain boundary conditions. 
In recent work by Gaspard et al., these methods have successfully been applied to the theory of dynamical systems 
pf|,p5 32 5?|,7^(|. In the following, the principles of first passage for the class V of dynamical systems defined above 
will be briefly outlined. The method will turn out to provide a convenient starting point for computing parameter- 
dependent diffusion coefficients. 

One may distinguish three different steps in applying the method: 

Step 1: Solve the one-dimensional phenomenological diffusion equation 

dn d 2 n , 

with suitable boundary conditions, where n :— n(x,t) stands for the macroscopic density of particles at point x and 
time t. This equation serves here as a definition for the diffusion coefficient D. 
Step 2: Solve the Frobenius- Perron equation 

p n +i{x) = dy p n (y) 8(x - M a (y)) , (8) 
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which represents the continuity equation for the probability density p n {x) of the dynamical system M a (y) |66j]8l|] . 

Step 3: For a chain of boxes of chainlength L, consider the limit chainlength L and time n to infinity: //for given 
slope a the respective largest eigenmodes of n and p turn out to be identical in an appropriate scaling limit, then 
D(a) can be computed by matching the eigenmodes of the probability density p to the particle density n: For periodic 
boundary conditions, i.e., n(0,t) = n{L,t) and p n (0) = p n {L), one obtains 

D(a) = (J^J ldec {a) , (9) 

where ^decio) is the decay rate in the closed system to be calculated directly from the Frobenius-Perron equation and 
therefore determined by quantities of the deterministic dynamical system. For absorbing boundary conditions, i.e., 
n(0, t) = n(L, t) = and p n (0) = p n {L) = 0, the same procedure leads to 

L" 2 



D(a) = lim - lesc {a) , (10) 
L^oo \ ir J 

where 7 esc (a) is the escape rate for the open system. This quantity can be further determined by the escape rate 
formalism P,p2fl to 

7 esc (a) = \(K;a) - h KS (U;a) , (11) 

where the Lyapunov exponent \(1Z; a) and the Kolmogorov-Sinai (KS) entropy hxsi'R-', a ) are defined on the repeller 
1Z of the dynamical system. This equation is an extension of Pesin's formula to open systems, which is obtained in case 
of 7 esc = 0. Eqs.(^) and ( |Io| ) have been applied to a variety of models, like the periodic Lorentz gas, two-dimensional 
multibaker maps and certain one-dimensional chains of maps, by Gaspard and coworkers |^,^] . Eq. ( |l0| ) , together 
with Eq. ( |TT|) , has first been presented for the two-dimensional periodic Lorentz gas |H] and has later been generalized 
to other transport coefficients and dynamical systems ]79|j80| ]. However, although of fundamental physical importance, 
it seems in general to be difficult to use this equation for practical evaluations of D(a), because usually the KS entropy 
is hard to calculate Jog ], Instead, Eq.(|lO|) with Eq.(pd|) can be inverted to get the KS entropy via the decay rate of 
the dynamical system of Eq.(||) to hKs{TZ;a) = \(lZ;a) — j^decifl) (L — > oo) , or by employing Eq. (||) via the 
diffusion coefficient in the limit of large L. 



III. SOLUTION OF THE FROBENIUS-PERRON EQUATION 



Following the first passage method, the problem of computing parameter-dependent diffusion coefficients essentially 
reduces to solving the Frobenius-Perron equation for the dynamical system in a certain limit. In this section, a general 
method will be presented by which this goal can be achieved. Its principles will be illustrated by performing analytical 
calculations for some special cases of parameters of map C Our method is based on finding Markov partitions and on 
defining respective transition matrices. This approach is quite well-known, especially in the mathematical literature 
p3| pot and has been employed by many authors for the calculation of dynamical systems quantities ]6^ , [70| , pl| -p9| . 
We apply it here for the first time to compute the full parameter-dependent deterministic diffusion coefficient. 



A. Transition matrix method 



As a first example, the diffusion coefficient D(a) shall be computed for map C at slope a = 4, as sketched in Fig. ||, 
supplemented by periodic boundary conditions. The calculation will be done according to the three-step procedure 
outlined above. 

Step 1: The one-dimensional diffusion equation Eq. (|7j) can be solved with periodic boundary conditions straight- 
forward to 

/ \ v ^ / .2-^771.5, \ / ,2-Km . 2i{m \ 

n{x,t) = a + ex P I -(— —) Dt I I a m cos(— —x) + b rn sm(— — x) I , (12) 

m=l \ / \ / 

where ao , a m and b m are the Fourier coefficients to be determined by an initial particle density n(x,0). 

Step 2: To solve the Frobenius-Perron equation, the key idea is to write this equation as a matrix equation p2| , |64[ . 
For this purpose, one needs to find a suitable partition of the map, i.e., a decomposition of the real line into a set of 
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subintervals, called elements, or parts of the partition. The single parts of the partition have to be such that they do 
not overlap except at boundary points, which are referred to as points of the partition, and that they cover the real 
line completely |64|]. In case of slope a = 4, such a partition is naturally provided by the box boundaries. The grid of 
dashed lines in Fig. |^ represents a two-dimensional image of the one-dimensional partition introduced above, which 
is generated by the application of the map. 

Now an initial density of points shall be considered which covers, e.g., the interval in the second box of Fig. || 
uniformly. By applying the map, one observes that points of this interval get mapped two-fold on the interval in 
the second box again, but that there is also escape from this box which covers the third and the first box intervals, 
respectively. Since map C is old, this mechanism applies to any box of the chain of chainlength L, modified only by 
the boundary conditions. Taking into account the stretching of the density by the slope a at each iteration, this leads 
to a matrix equation of 



Pn+l 



-T(a)p„ 
a 



(13) 



where for a — 4 the L x i-transition matrix T(4) can be constructed to 



T(4) 



/2 1 ••• 1\ 
12 1 ••• 
1 2 1 ••• 



••• 1 2 10 
••• 1 2 1 
VlO ••• 21/ 



(14) 



The matrix elements in the upper right and lower left edges are due to periodic boundary conditions and reflect the 
motion of points from the Lth box of the chain to the first one and vice versa. 

In Eq. dig)) , the transition matrix T{a) is applied to a column vector p n of the probability density p n (x) which, in 
case of a = 4, can be written as 



Pn = \pn{x) >:= (pl,pl 



.Pn.- 



(15) 



where "*" denotes the transpose, and p„ represents the component of the probability density in the fcth box, p n (x) — 
p\ , k — 1 < x < k , fc = 1, p„ being constant on each part of the partition. 

In case of a = 4, the transition matrix is symmetric and can be diagonalized by spectral decomposition. Solving 
the eigenvalue problem 



T(4) \(j) m {x) >= Xm (4) \4> m {x) > , 
where Xm(4) and \<j) m {x) > are the eigenvalues and eigenvectors of T"(4), respectively, one obtains 



(16) 



L-l 



\p n (x) > = - ^ Xm(4) \4> m {x) >< (f> m (x)\p n (x) 



> 



m=0 
L-l 



exp ( —n In 



Xm(4) 



\4> m (x) X (j) m {x)\p (x) > 



(17) 



where IpoO^) > is an initial probability density vector and In 4 is the Lyapunov exponent of the map. Note that the 
choice of initial probability densities is restricted by this method to functions which can be written in the vecto r form 
of Eq.(|l5|). For matri ces of the type of T(4), it is well-known how to solve their eigenvalue problems [ 100| -102|. This 
is performed in App. A 1 for a more general case, which includes the example under consideration. For slope a — A, 
one gets 



X m (4) = 2 + 2cos6> r , 
\<f> m (x) > = (^,$ 



A- 



An 



27T 

— m , m — 0, . 



,4>m = a>m<t> k m , 1 + b m <f>t 2 



cos 9 m (k ~ 1) , </^ 2 := sin6> m (fc - 1) 
k = 1 , . . . , L , k — 1 < x < k 



(18) 
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with a m and b m to be fixed by suitable normalization conditions. 

Step 3: To compute the diffusion coefficient D(4), it remains to match the first few largest eigenmodes of the 
diffusion equation to the ones of the Frobenius-Perron equation: In the limit as the time t and the system size L 
approach infinity, the particle density n(x,t), Eq. (|l2|) , in the diffusion equation becomes 



n{x, t) ~ const. + exp 



2tt 2tt n 

Acos( — x) + B sm(—x j 
L L 



(19) 



where the constant represents the uniform equilibrium density of the equation. 

Analogously, for discrete time n and chainlength L to infinity, one obtains for the probability density p n {x) of the 
Frobenius-Perron equation, Eq.(^) with Eq.(|i8|), 



2ir ~ 2ir 

p n (x) ~ const. + exp (-7d ec (4)n) ( Acos(— (k - 1)) + Bsin(— (k - 1)) 



k = 1, . 



with a decay rate of 



k - 1 < x < k 



7dec(4) = In- 



2cos^ 



(20) 



(21) 



of the dynamical system, determined by the second largest eigenvalue of the matrix T(4), see Eq.(|l8|). Note that the 
largest eigenvalue is equal to the slope of the map so that for the first term in Eq.(|2^) the exponential vanishes, and 
one obtains a uniform equilibrium density. Apart from generic discretization effects in the time and position variables, 
which may be neglected in the limit of time to infinity and after a suitable spatial coarse graining, the eigenmodes of 
Eqs.([l9|) and (p0|) match precisely so that, according to Eq.(0), the diffusion coefficient -D(4) can be computed to 



D(4) 



(^) 2 7dec(4) = 



0(L~ 



This result is identical to what is obtained from a simple random walk model |4qj58|j7£ 
generalized straightforward to all even integers values of the slope, as is shown in App ^A 1\ , 
dependent diffusion coefficient of 



D(a) = ^(a-l)(a-2) 



= 2k 



k e n 



(22) 

j. The procedure can be 
and leads to a parameter- 



(23) 



in agreement with results of Ref. E[ . 

A slightly more complicated example is the case of slope a — 3, see, e.g., Fig. |2[ which will be treated in the 
following. Analogously to the previous example, for a=3a simple partition can be constructed, the parts of which 
are all of length 1/2. According to this partition, a transition matrix T(3) can be determined, given schematically by 



T(3) 





1 
















1 


\ 


1 


1 





1 
















1 





1 


1 






















1 


1 





1 
















1 





1 


1 










\o 


1 
















1 


1 / 



(24) 



Note that, in contrast to the case of a — 4, here the matrix is formed by submatrix blocks which move periodically 
to the right every two rows. Since the partition of a = 3 is a bit more complicated than for a = 4, the blocks refer 
to the partition of each box, whereas the shift again is related to the lift property of the old map. The matrix T(3) 
is not sy mmet ric. However, the eigenvalue problem of this matrix can still be solved analogously to the case of a = 4 
(see App. AT). The s pectr um of the matrix turns out to be highly degenerate, and therefore T(3) cannot be simply 
diagonalized anymore [103]. This is due to the fact that the matrix T(3) is non-normal, i.e., T(3)T*(3) ^ T*(3)T(3), 
which means that it does not provide a system of orthogonal eigenvectors. Of course still a transformation onto Jordan 
normal form to "block-diagonalize" this matrix could be applied. However, thi s see ms to be more useful for proving 
mathematical theorems than for analytical calculations of diffusion coefficients [104]. Since the probability density of 
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the Frobenius-Perron equation is determined by iteration of transition matrices, cf. Eq. (|l3|), for slope a = 3 it is not 
in advance clear how the single eigenmodes "mix" in the limit time and chainlength to infinity and whether they can 
be matched to the solutions of the diffusion equation as before so that the diffusion coefficient D(3) is again simply 
determined by the second largest eigenvalue of the matrix. 

This problem will first be approached pragmatically. Analogously to the analytical solutions of Eqs. (|f^) and ( pc| ) 
for slope a = 4, Fig. [| shows a plot of the two second largest eigenmodes of T(3) in comparison to the solution of the 
diffusion equation. Again, one observes total agreement, except differences in the fine structure. The same is true 
for the other first few largest eigenmodes of T(3). Thus, although straightforward diagonalization and, therefore, a 
simple solution of the Frobenius-Perron equation like Eq. ( |l7| ) are not possible anymore, the largest eigenmodes of 
T(3) behave correctly in the sense of the phenomenological diffusion equation so that it is suggestive to compute the 
diffusion coefficient D(3) via the second largest eigenvalue of T(3) again. With 



see App. Al, and Eq.(^), one gets 



7(3) = In 



1 + 2 cos(2tt/L) ' 



(25) 



D(3) = - + 0(L-*) . 



(26) 



As for a = 4, this result is obtained as well from a simple random walk model. However, to produce this value, the 
respective random walk has to be defined in a slightly different way than for a — 4 58 7q]. Analogously to the case 
of ev en in teger slopes, the exact calculations can be generalized to all odd integer values of the slope and lead to (see 
App. |AT| ) 



= 2k-l 



k e n 



(27) 



which again is identical to the result of Ref . . 

Before this approach will be extended to other parameter values of the slope in the following section, some details 
will be discussed which are closely related to the calculations above: A special feature of the diffusion coefficient 
results for integer slopes shall be pointed out; further limits of this method with respect to generalized diffusion 
coefficients shall be critically discussed; and the application of the method to absorbing boundary conditions shall be 
briefly outlined. 



1. Diffusion coefficients for integer slopes 



Eqs. (|22j ) and (|26j) show already that _D(4) < -D(3), which is at first sight counterintuitive. By evaluating the general 
formulas of D(a) given by Eqs. ( |23] ) and (|27]) at other even and odd integer slopes, one realizes that this inequality 
reflects a general oscillatory behavior of D[a) at integer slopes. This result has already been obtained by Fujisaka 
and Grossmann [ffl|| , but it has not been discussed in further detail in their work. A sim i lar o scillatory behavior has 
been observed for deterministic diffusion in certain classes of two-dimensional maps 105 - IPS) ] . This behavior cannot 
be understood completely by one consistent simple random walk model |5qf7q|. 



2. Matching lower eigenmodes 

There appear serious problems in trying to extend the matching eigenmodes procedure to arbitrarily low eigenmodes, 
even in case of a = 4, where the matrix is diagonalizable. With Eq. (|l8|) , one can check that 

<t>m,l = 4>L-m,l : 4L,2 = -4> h L ~ m ,2 i k = l,...L , m = 1, ...,£- 1 , (28) 

i.e., in contrast to the m eigenmodes of the diffusion equation Eq. ([?]), the frequency of the eigenmodes of T(4) does 
not increase monotonically with m, but gets "flipped over" at the (£/2)th (for L even) eigenmode such that the first 
and the last half of the number of eigenmodes are identical, except a minus sign. This is due to the discretization of 
the position variable x in the diffusion equation to k in the Frobenius-Perron matrix equation Eq. (|l^) , which was one 
of the basic ingredients for the possibility to construct transition matrices. 

Moreover, one should note that, according to Eq.([l8|), the sm alles t eigenval ue o f T(4) is equal to zero. For T(3), 
a large number of eigenvalues are even less than zero (see Eq. ( A16| ) in App. A 1 for a = 3). Thus, except for the 
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first few largest eigenmodes, which still match reasonably well to the eigenmodes of the diffusion equation in the limit 
time n and chainlength L to infinity, one cannot expect the method to work simply that the components of a "time- 
dependent" diffusion coefficient D n {a) are determined by smaller eigenvalues of the transition matrices in straight 
analogy to Eq.(||). This could be taken as a hint that, to obtain more details of the dynamics, refined methods are 
needed. For example, in Ref. ]ll(fl the first orders of a position-dependent diffusion coefficient have been determined 
for a class P-map according to a procedure which avoids the discretization of the real line. 



3. Absorbing and periodic boundary conditions 



The same procedure as outlined for periodic boundary conditions can also be employed for absorbing boundaries. 
It shall be sketched briefly, according to the three steps distinguished before: 

Step 1 : The one-dimensional diffusion equation with absorbing boundary conditions can be solved to 

n(x, t) = ^a m exp (-( — ) Dtj sm(—x) (29) 

TCI — 1 

with a m denoting again the Fourier coefficients. 

Step 2 : The transition matrices for a = 4 and a = 3 at these boundary conditions are identical to the ones of 
Eqs.([l4|j",(|2~4|), except that the matrices now contain zeroes as matrix elements in the upper right and lower left 
corners. However, due to this slight change in their basic structure there is no general method to solve the eigenvalue 
problems for this type of matrices anymore, in contrast to the case of periodic boundary conditions. At least for 
a = 3 and a = 4, it is still possible to obtain analytical solutions by straightforward calculations analogous to the 



ones performed in Ref. |32j (see App. A 2), but for any higher integer value of the slope even these basic methods fail. 
This appears to be caused by strong boundary layers. Fig. ^| shows numerical solutions for the largest eigenmodes of 
the first odd integer slope transition matrices in comparison to the solutions of the diffusion equation Eq. (fjj) (details 
of the numerics applied here are given in the following section). It can be seen that near the boundaries, there 
are pronounced deviations between the Frobenius-Perron and the diffusion equation solutions. These deviations are 
getting smaller in the interior region of the chain, but are gradually getting stronger with increasing the value of 
the slope, as is shown in the magnification. The same behavior can be found for even integer slopes, although the 
quantitative deviation of these eigenmodes from the ones of the diffusion equation solutions is slightly less than for odd 
values of the slope. Thus, obviously absorbing boundary conditions disturb the deterministic dynamics significantly, 
whereas similar effects do not occur for periodic boundary conditions, which therefore could be characterized as a 
kind of "natural boundary conditions" for this periodic dynamical system. 

Step 3: The different boundary conditions do not only show up in the eigenmodes of the transition matrices, but 
also in the calculation of the diffusion coefficients. In analogy to periodic boundary conditions, the escape rate of the 
dynamical system at a = 3 and a = 4 is determined to 

7 esc (a) = In — — (30) 

Xmax [&) 



with Xmax (a) being the largest eigenvalue of the transition matrix (see App. A 2) 



; (3) = 1 + 2 cos and Xmax (4) = 2 + 2 cos . (31) 



Feeding this into Eq.([l0|) via matching eigenmodes, one obtains 



which gives a convergence of the diffusion coefficient with the chainlength L significantly below that obtained from 
periodic boundary conditions. For example, for a chainlength of L = 100 the convergence is about two orders of 
magnitude worse. It can be concluded that the transition matrix method works in principle for absorbing bound- 
ary conditions as well, but that here its range of application to compute diffusion coefficients is qualitatively and 
quantitatively more restricted because of long-range boundary layers. 
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B. Markov partitions 



In the previous section, the choice of simple partitions enabled the construction of transition matrices. These 
matrices provided a way to solve the Frobenius-Perron equation in a certain limit. However, so far this method has 
only been applied to very special cases of map £, defined by integer slopes. This raises the question whether an 
extension of this method to other values of the slope is possible. For this purpose, the idea of choosing a suitable 
partition of the map has to be generalized. Taking a look at Fig. |§| again, one observes that the graph of the map 
"crosses" or "touches" a vertical line of the grid only at some grid points. Furthermore, the local extrema of the map, 
which are here identical to the points of discontinuity, are situated on, or just "touch" horizontal lines of the grid, 
whereas other crossovers of horizontal lines occur at no specific point. The same characteristics can be verified, e.g., 
for the respective partition of slope a = 3. These conditions ensure that it is possible to obtain a correct transition 
matrix from a partition, since to be modeled by a matrix, a density of points, which covers parts of the partition 
completely, has to get mapped in a way that its image again covers parts of the partition completely, and not partially. 
This basic property of a "suitable partition to construct transition matrices" is already the essence of what is known 
as a Markov partition: 

Definition 1 (Markov partition, verbal definition) For one- dimensional maps, a partition is a Markov partition 
i f and only if parts of the partition get mapped again onto parts of the partition, or onto unions of parts of the partition 

A mo re fo rmal definition of one-dimensional Markov partitions, as well as further details, can be found in Refs. 
|p7,30 1 11 1 . The next goal must be to find a general rule of how to construct Markov partitions for map £ at other, 



non-trivial parameter values of the slope. Because of the periodicity of the chain of maps it suffices to find a Markov 
partition for a single box map, i.e., for the respective map in one box without applying the modulus to restrict it onto 
the unit interval. Here, the fact can be used that the extrema, which are the critical points of the box map, have to 
touch horizontal lines, as explained before, which means that to obtain a Markov partition the extrema have to get 
mapped onto partition points. Since any box map of map C is symmetric with respect to the point (x, y) = (0.5, 0.5), 
the problem reduces to considering only one of the extrema in the following, e.g., the maximum. Changing the height 
of the maximum corresponds to changing the slope of the map. Therefore, if one needs to find Markov partitions for 
parameter values of the slope, one can do it the other way around by the following Markov condition: 

Definition 2 (Markov condition, verbal definition) For map C, Markov partition values of the slope are de- 
termined by choosing the slope such that the maximum of the box map gets mapped onto a point of the partition 
again. 

In Fig. |^, four examples of non-trivial Markov partitions for map C are depicted with respect to their box maps. One 
may check that the handwaving conditions which have been extracted from the integer slope examples to motivate 
Definition [l] are fulfilled and especially that the partitions shown in the figure obey the Markov partition definition 
|l|. The detailed structure of the partitions can be arbitrarily complex. In Fig. || (a) and (b) a special orbit has been 
marked by bold black lines with arrows. It represent what will be called the generating orbit of a Markov partition: 
For the two examples shown here, the starting point of this orbit is given by the maximum of the preceding box map, 
since this maximum must also be a partition point. The iterations of this orbit, as indicated by the arrows in the 
figure, define the single partition points. This way, the number of partition parts is related to the number of iterations 
of the generating orbit. In case of Fig. [j| (a) and (b) the orbit is eventually periodic, i.e., it finally gets mapped onto 
the fixed point x — 0, however, it can also be periodic with a certain period, as, for example, in case of Fig. ^ (d) 
(period four). Thus, the generating orbit is the key of how to find Markov partition values of the slope in a systematic 
way. 

On this basis, a general algebraic procedure to compute such values of the slope can be developed. One starts with 



a further topological reduction of the whole chain of boxes |112|. Since map C is old, it is possible to construct the 
Markov partition for the whole chain from a reduced map 

M a (x) := M a {x) mod 1 (33) 

via periodic continuation, where x := x — [x] is the fractional part of x, x £ (0, 1], and [x] denotes the largest integer 
less than x. Therefore, it remains to find Markov partitions for map M a (x) of the equation above. This can be done 
in the following way: Let 

M a (i),l-M a (i)| , e<~ , (34) 







be the minimal distance of a maximum of the box map M a (x) to an integer value. With respect to the Markov 
condition given by Definition ^], it is clear that e has to be a partition point. Since M a (x) is point symmetric, 1 — e 
also has to be a partition point, and because of map C being old, the fixed point x — is necessarily another partition 
point. Thus, the reduced map governs its internal box dynamics according to 

x n+ i = M a (x n ) , x n = M"(x) , x = x . (35) 

Since 0, e and 1 — e have to be partition points, the Markov condition Definition || can be formalized to 

M a "(e)=<5 , 5e{Q,e,l-e} , (36) 

i.e., the generating orbit of a Markov partition is defined by the initial condition e, its end point S and the iteration 
number n. According to Eq.(|34|), e is determined by the slope a. Therefore, for map C Markov partition values of the 
slope can be computed as solutions of Eq.(|36|). The evaluation of this equation can be performed numerically as well 
as, to a certain degree, analytically. To obtain analytical results for Markov partition values of the slope, one has to 
determine the structure of the generating orbit in advance, i.e., one has to know whether it hits the left or the right 
branch of the box map at the next iteration. Then one can write down an algebraic equation which remains to be 
solved. For example, for the Markov partition Fig. ^ (c), the generating orbit is determined by 

x 1 = M a (e) , e<i and 6 = M a { Xl ) - 3 , x x < ^ (37) 

with a = 2(3 + e) and 5 = at iteration number n = 2. This leads to 

a 3 - 6a 2 - 6 = , a > 2 , (38) 

for which one may verify a ~ 6.158 as the correct solution. This way, all Markov partition values of the slope are 
the roots of algebraic equations of (n + l)th order. More examples with analytical solutions are discussed in App. 



A3. Since one usually faces the problem to solve algebraic equations of order greater than 3, numerical solutions of 
Eq.(|36"|) are desirable, although one should take into account that iterations of the reduced map M a (x) contain many 
discontinui ties, due to the original discontinuity of M a (x) at x = 1/2 as well as due to applying the modulus to M a (x) 



in Eq.(||) (TT| 



With respect to the three different end points S of the generating orbit in the formal Markov condition Eq.(p6[), 
three series of Markov partitions can be distinguished. For each series one can increase the iteration number n, and 
one can vary the range of the slope a systematically. These three series have been used as the basis for numerical 
calculations of the diffusion coefficient D(a), as will be explained below. However, there exist additional suitable end 
points S for the generating orbit. As an example, one can choose 5 to be a point on a two-periodic orbit, 

MZ(S)=5 , o = 2(l + c) , 0<e<I 5 = — i±^_ ( 39 ) 

so that the generating orbit is again eventually periodic, but now being mapped on a periodic orbit which is part of 
the Markov partition instead of being mapped on a simple fixed point. This way, certain periodic orbits can serve for 
defining an arbitrary number of new Markov partition series with respect to the choice of respective new end points 
S. On the other hand, the set of Markov partition generating orbits is not equal to the set of all periodic orbits. For 
example, for the range 2 < a < 3, Eq.([59|) shows that there exists a two-periodic orbit for any slope a, but not any 
maximum of the map in this range necessarily maps onto this periodic orbit, as is already illustrated by Fig. [| (a) 
and (b) , or by other simple solutions of Eq. (^6|) , respectively. This proves that Markov partition generating orbits are 
in fact a subset of all periodic orbits of the map. 

With respect to varying the iteration number n and the end point 5, one can expect to get an infinite number of 
Markov partition values of the slope. In fact, for certain classes of maps the existence of Markov partitions can be 
considered as a natural property of the map [^3|j85|,|8^,|90| . According to the explanations above, this does not seem 



to be true for map C. Instead, there is numerical evidence for the following conjecture [114| 



Conjecture 1 (Denseness property of Markov partitions) For map C, the Markov partition values of the slope 
a are dense on the real line with a > 2. 

This denseness conjecture should ensure that it is possible to obtain a representative curve for the parameter-dependent 
diffusion coefficient D(a) solely by computing diffusion coefficients at Markov partition values of the slope. Conjecture 
|l| may hold for all other class P-maps as well. To do such computations, one needs to construct the corresponding 
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transition matrices to the Markov partitions obtained, as it has been shown for the slopes a = 3 and a — 4 of map 
C. This can be done according to the following rule: Take as an example any of the box map Markov partitions 
illustrated in Fig. [| e.g., case (a). Any dashed rectangle of this partition may be denoted as a cell of the partition. 
These single cells correspond to the single entries, or matrix elements, of the transition matrix to be obtained. The 
transition matrix corresponding to this Markov partition can now be constructed by checking where the graph of the 
map goes across a cell of the partition, by counting the number of these occurrences in each cell and by writing down 
these values as the matrix elements. For map C, usually these matrix elements will consist of zeroes and ones, but 
the way they are defined here they can also take other integer values, depending on the choice of the partition, as, 
e.g., illustrated in case of a = 4, Eq. (p"l|). 

The construction of the box map transition matrix can be simplified by taking the point symmetry of the box map 
into account. The transition matrix of the full chain of chainlcngth L again follows by periodic continuation. These 
matrices can be denoted as topological transition matrices, since they reflect purely the topology of the map with 
respect to the Markov partitions, without involving any transition probabilities at this point. In Refs. | TO 8^j85| , p6| , p0[ | 
mathematically rigorous definitions of these transition matrices can be found. The property of map L being old 
induces a certain structure in the topological transition matrices. They are said to be banded square block Toeplitz 
matrices, i.e., they consist of certain submatrices, called blocks, corresponding to the box map Markov partitions, 
and these bloc ks are the s ame along diagonals of the topological transition matrix parallel to the main diagonal, 
forming bands [102.115 116|. Applying periodic boundary conditions to the chain of boxes defines a subclass of these 
Toeplitz matrices, called block circulants, where each row is constructed by cycling the previous row forward one 
block [ 100] - 102, 115 1, see, e.g., the matrices T(4), Eq.(14), as an example for a simple circulant and T(3), Eq.(p4[), 
for a block circulant. According to the transition matrix method outlined in the previous section, it remains to 
solve the eigenvalue problems of these matrices and to match the respective eigenmodes to the ones of the diffusion 
equation for computing the corresponding diffusion coefficients D(a). Here, periodic boundary conditions are of great 
advantage. Analytically, as mentioned befor e and as shown in App. A 1, there exists a general procedure how to solve 
the eigenvalue problems of simple circulants 1 100 102 1 , and in s ome case s it is possible to reduce the eigenvalue problem 



of a block circulant to that of a simple circulant (see App. Al and A3). If this method works, it automatically yields 
"nice eigenmodes", i.e., eigenvectors of the form of sines and cosines with so me fi ne structure. These eigenmodes are 
similar to the eigenmodes of the diffusion equation at this stage (see App. A3), i.e., before iterating the matrices 
according to the Frobenius-Perron matrix equation Eq. (|l3|) . The situation is quite different for absorbing boundary 
conditions, where no such general procedure exists (see App. A 2). 

If analytical solutions of the eigenvalue problems are not possible anymore, one can obtain numerical solutions. 
Well-known software packages like NAG and IMSL provide subroutines to solve the eigenvalue problems of these 
matrices. Unfortunately, the numerically obtained results for the full spectra turned o ut n ot to be very reliable to 
a certain extent: In comparison to analytical results for periodic boundaries (see App. A3), the NAG package does 
not compute all eigenvectors correctly, i.e., in the numerical results usually some linear independent eigenvectors are 
missing. Moreover, both packages provide spectra of eigenvalues which, although partly identical to the analytical 
solutions , dif fer in their full range quantitatively to the ones calculated analytically, not taking any degeneracy into 
117|. Such numerical problems seem to be inherent to the class of non-normal Toeplitz matrices, as has 

In fact, solving eige nvalue p roblems for Toeplitz matrices can 

However, solely for the purpose 
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account 

already been pointed out by Beam and Warming [ 1 1 5 1 . 
be considered as a field of active recent research in numerical mathematics 
of computing diffusion coefficients not the full spectra of the transition matrices are required, but only the few largest 
eigenvalues and eigenvectors are of interest. With respect to eigenvectors, the IMSL package has been checked to be 
reliable in this range, and with respect to eigenvalues, both packages provide exact and identical numerical results, 
especially for the second largest eigenvalue, which determines the diffusion coefficients. From a practical point of 
view, the NAG package is considerably more efficient in the use of computing time than IMSL and more flexible 
in its subroutines. Therefore, for quantitative computations of eigenvalues and diffusion coefficients NAG routines 
have been used, whereas for eigenvector computations IMSL procedures have been employed. For computations 
of diffusion coefficients, it is also favorable to consider only the case of periodic boundary co nditio ns, i.e., solving 
eigenvalue problems for block circulants, respectively, since it has already been discussed in Sect III A that absorbing 
boundaries lead to a poor convergence rate of the diffusion coefficient with the chainlength L. 

Fig. ^ contains two examples of second largest eigenmodes for chains of boxes with periodic boundaries and non- 
trivial Markov partitions. Again, one gets "nice" second largest eigenmodes, i.e., functions which behave like sines 
and cosines on a large scale. However, the structure of these eigenmodes is much more complex on a fine scale, as 
can be seen in the magnifications of certain regions. The periodic continuation of the fine structure suggests that it is 
related to the dynamics of the box map, and therefore varies with changing the slope, whereas the general large-scale 
behavior of the eigenmodes seems to be a property of the chain of boxes which shows up independently from such 
microscopic details. These characteristics have been checked numerically for a variety of other Markov partition 
values of the slope and seem to be a universal feature of map L, and probably of all class "P-maps. One may assume 
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that the fine structure is somehow related to the strength of the diffusion coefficient and that, on the other hand, 
the universal large-scale structure of the eigenmodes is related to the existence of diffusion coefficients for non-trivial 
Markov partition values of the slope. In fact, the specific character of the eigenmodes discussed above, which shows 
up in any analytical solution of (block) circulants and which is supported by numerical results, forms the basis for 
the following conjecture: 

Conjecture 2 (Existence of diffusion coefficients) Let M a (x) be a class V-map. If for given value of the slope 
the map is uniquely ergodic and if there exists a Markov partition, then the map is diffusive. 

To our knowledge, so far no proof has been given in the literature for the existence of diffusion coefficients in class 
"P-maps for general value of the slope. However, dealing with a rigorous foundation of the transition matrix method 
turns out to be intimately connected to proving the existence of diffusion coefficients in this class of dynamical 
systems. Without going into too much detail here, some remarks are in order to provide at least a motivation for 
this conjecture: The existence of Markov partitions guarantees that exact transition matrices can be used. The 
restriction to class "P-maps ensures that topological transition matrices can be constructed in the simple way outlined 
before, and the old property included in the definition of class V determines the global structure of the topological 
transition matrices such that the eigenmodes are "nice", at least for periodic boundary conditions. The requirement 
to be uniquely ergodic establishes the possibility of diffusion in the chain of boxes and confirms also the uniqueness 
of the diffusion coefficient to be obtained (a simple counterexample shows that not any chain of boxes with escape 
out of one box is automatically diffusive) . Finally, the term diffusive shall be understood in the sense that a diffusion 
coefficient exists as defined by the statistical diffusion equation Eq.(Q), which has been introduced to the dynamical 
system by successfully performing the matching eigenmodes procedure outlined in the previous section. Therefore, the 
main proposition of this conjecture is that the matching eigenmodes procedure required by the first passage method 
works for any value of the slope, if the respective conditions are fulfilled. A corollary to this conjecture is that in 
the limit of time n and chainlength L to infinity, the Frobenius-Perron equation of the respective class P-dynamical 
systems always provides "nice", i.e., the correct diffusive eigenmodes. As another corollary, it follows that there is no 
anomalous diffusion in class P-maps, i.e., that normal diffusion is "typical" for such piecewise linear maps. A rigorous 
math ematical proof of Conjecture g seems to be possible along the lines of first passage and transition matrix method 



[104] 



Results base d on this method shall be presented in the next section. They have been verified by an other numerical 
method |5^ , |77| , |l04"| ] , by another analytical method which ha s be en implemented numerically [58 120| , as well as, to a 



certain degree, by straightforward com pute r simulations |38 104]. Meanwhile, the same results have been obtained by 



J. Groeneveld with a different method [121], and they have also partly been reproduced by cycle expansion techniques 



IV. FRACTAL DIFFUSION COEFFICIENTS: RESULTS 



Based on the methods presented in the previous section, the parameter-dependent diffusion coefficient D(a) has 
been computed numerically for map C for a broad range of values of the slope. The main results are shown in Fig. 
||. The numerical precision obtained depends on the convergence of the diffusion coefficient with the chainlength, 
cf. Eq.(32), and is better than 10 -4 for each D(a) so that error bars do not appear in the diagrams. It should 
be emphasized that the numerical method employed here was the first one by which these curves of D(a) have 
been obtained. It is by far not the mos t efficient one of the procedures developed up to now to compute deterministic 
diffusion coefficients 
transport coefficients 



77 104 12C . 121 ] . However, it turns out to be very useful as well to compute other deterministic 

where more efficient methods fail. 



_., chemical reaction rates |3S 
Fig. [§] (a) shows the diffusion coefficient of map C for values of the slope in the range 2 < a < 8. The strength of 
diffusion clearly increases globally by increasing the slope from a = 2 to a = 8. This might be expected intuitively, 
since the probability of a particle to escape out of a box, as well as the mean distance a particle travels by performing 
a jump, are getting larger by increasing the value of the slope J|j|(78|. However, the increase of the diffusion coefficient 
is not monotonic and consists of oscillations not only at integer values of the slope, as has already been mentioned 
in Section EL but also on much finer scales between integer values. In fact, Fig. | (a) shows a certain regularity in 
the appearance of "wiggles", i.e., local maxima and minima. If one denotes the local maxima at odd integer slopes 
a as wiggles of 0th order and any smaller local maxima systematically as wiggles of higher order, one can find one 
maximum of first order below a = 3, three maxima of first order in the range 3 < a < 5, five maxima of first order 
in the range 5 < a < 7, .... This regularity even persists to a certain extent on finer scales, although according to 
a slightly different rule, as can be seen, e.g., in the magnification Fig. || (f), 6 < a < 7, where exactly six wiggles of 
second order appear between the repsective wiggles of first order. In the same way, six wiggles of third order can be 
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observed in this region in further magnifications, and similar structures show up in the region of 4 < a < 5 with four 
wiggles of second and four wiggles of third order. The region of 2 < a < 3 is somewhat special and will be discussed 
separately. Thus, while the number of wiggles of first order increases by a step of two with increasing the slope, the 
number of wiggles of higher order remains constant in the region between two respective wiggles of first order, even by 
increasing the order of the wiggles to be considered. On the other hand, magnifications of other regions of the slope 
show that the structure of the curve is not everywhere that simple. For example, blow-ups of the regions 3 < a < 4, 
Fig. H (b), and 5 < a < 6, Fig. || (d), do not enable a clear distinction between "wiggles of different orders" anymore. 
Instead, they provide more complex structures which further magnifications, as, e.g., Figs. |§| (c) and (e), reveal to be 
self-similar. 

It can be summarized at this stage that different regions of the curve exhibit different kinds of self-similarity, partly 
being fairly simple, but partly also being highly non-trivi al. T hus, the results of Fig. |^ suggest that the parameter- 
dependent diffusion coefficient D(a) for map C is fractal 122 1. More evidence for the fractality of the curve can be 
obtained in three different ways: Firstly, qualitative and quantitative explanations for the wiggles in certain regions 
of the slope will be provided, which ensure that these regions exhibit non-trivial self-similar behavior. This will be 
demonstrated in the following. Secondly, it is striking to observe that esp ec ially diagrams (c), (e), and (f) resemble 
graphs of some fractal functions, which have been obtained in Refs. |33||34| , |l23| ] by working on dynamical systems 
being v ery s imilar to the ones considered here. These functions have been shown to possess fractal dimensions close 
to one [123 1 . On the one hand, this gives further evidence for the fractality of the curves of Fig. ^. On the other 
hand, this raises the new question whether there is an analytical fractal representation for certain regions of the curve, 
or maybe e ven f or the full parameter-dependent diffusion coefficient of map C This problem is discussed in detail 
in Refs. [ |58|l20| . Thirdly, it should be mentioned that numerical computations of the box counting dimension (or 
capacity) ]66|] of the curve have been performed. The results indicate that the curves shown in Fig. |8| (a) - (f) have 
fractal dimensions d very close to, but not equal to one in a range of d = 1 + A<i , < Ad < 10 -2 . Because of the 
limited data set better values are difficult to get, especially since the fractal dimension is expected to be close to one 
in this case. Nevertheless, more detailed investigations of the fractal dimension for various regions of the curve could 
be of much interest. With respect to the magnifications in Fig. |[ it may even be conjectured that the full £>(a)-curve 
has locally different values of fractal dimensions. 

Fig. H illustrates the principles of a first qualitative approach to understand the occurrence of wiggles of Oth and 
1st order. It will be called plus-minus approach. The basic idea of this approach is to establish a connection between 
the appearance of wiggles in the -D(a)-graphs and the occurrence of certain dynamical correlations in the chain of 
boxes. These correlations are a main feature of transport of particles from one box to another, and they show up 
and disappear by varying the slope of the map. In the following, particles will be referred to solely by their positions, 
i.e., by points on the real line. Fig. ^| (a) and (b) sketch correlations of Oth order: As a starting point, the escape 
of particles out of one box in one direction, i.e., to the right, will be considered for varying the slope in the range 
2 < a < 4. Such an escape of points is related to a certain subinterval of the box which will be called escape region, 
as is shown in the figure. If points get mapped to the right at the next iteration, the respective subinterval will be 
denoted with a plus sign. The same way, subintervals will be denoted with a minus if points get mapped to the left. 
Therefore, the escape region marked in Fig. ^| is part of a plus region, and for small enough slope after only one 
iteration points of it get mapped directly into another plus region. This enhances diffusion, since particles can move 
continuously in one direction, i.e., here to the right. The behavior persists for increasing the slope up to a = 3. For 
slopes above this value, an increasing number of points of the escape region is now mapped into the minus region of 
the next box. This way, one obtains a "plus-minus" correlation, which means that particles either get slowed down, or 
even get scattered back into the previous box at the next iteration, which is surely bad for obtaining a strong diffusion 
coefficient. This game can be played by gradually increasing the value of the slope and leads to the qualitative "curve" 
in Fig. ^| (b) , which explains the oscillations at integer slopes and the wiggles of Oth order, respectively. The sequences 
which mark the extrema in this graph give the symbolic dynamics of orbits close to, but less than x — 1/2 after one 
iteration with respect to the reduced map Eq.(p3|), where the region < x < 1/2 has been labeled with a plus and 
1/2 < x < 1 with a minus. 

In Fig. H (c) the number of iterations has been increased to two. The method is the same as explained before, 
however, a further distinction has been made after the first iteration: new subintervals have been defined, which 
refer to points of the escape region being mapped to another plus or minus region at the second iteration. One can 
see that increasing the slope corresponds to creating different plus-minus sequences for orbits close to, but less than 
x = 1/2. This leads to the particles being in a good or bad position for going further in one direction with respect 
to the next iteration, depending on the value of the slope. The -D(a)-graph in Fig. || (d) again gives the qualitative 
behavior of D(a) to be expected with respect to the dynamical correlations after two iterations, up to a = 5. This 
result corresponds well to the number of wiggles of first order estimated in the respective regions of the slope. Again, 
the plus-minus sequences give the symbolic dynamics of points close to, but less than x — 1/2 after two iterations. 

The plus-minus method works on this level as well for any higher values of the slope and leads to a qualitative 
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explanation of the number of wiggles of first order for any region of the slope. To a certain degree, it can even 
explain additional features of the structure of the Z?(a)-curves: For example, in Fig. || (b) one observes that the local 
maximum is not precisely at a = 3, although this could be expected from the results of the plus-minus method of Oth 
order. Actually, particles of the escape region close to x = 1/2 can still reach a good position for further movement 
in one direction, even for slopes slightly above a = 3. This is due to the fact that, although such points get first 
scattered back into the previous box after two iterations, here they are now in an excellent position for further jumps 
to the right again. This way, these orbits perform a kind of "spiral" and seem to be responsible for the surprising fact 
that the odd integer slope values of D(a) are not precisely identical to the local extrema of the curve, but that there 
is always a kind of overhang, i.e., a further increase of the diffusion coefficient right above odd integer slopes, as, e.g., 
is shown in detail in Fig. |s| (c). 

Although the plus-minus method can be applied to achieve a qualitative understanding of the wiggles of Oth and 
1st order, further refinements of this method to obtain wiggles of higher order generally turned out not to be very 
promising. The main reason is that in case of more iterations of points of the escape region, the dynamics is getting 
quite complicated and is not easy to handle anymore in the qualitative way illustrated in Fig. ^|. However, the basic 
idea of this method can be made more quantitative by a procedure which shall be called turnstile dynamics. The 
principle of turnstile dynamics is again to investigate the appearance and disappearance of long-range dynamical 
correlations by iterating points with respect to varying the slope. The new feature is now that not the full intervals 
of all single boxes are taken into account. Instead, the analysis is restricted solely to the regions of the boxes where 
transport of particles from one box to another occurs in form of jumps. These regions are called turnstiles: 

Definition 3 (turnstile) Turnstiles are the "coupling regions" of the single boxes of a chain of class V, where points 
of one unit interval get mapped outside that particular interval into another unit interval. 

This notation has been adapted from the theory of transport of two-dimensional twist maps, such as sawtooth maps, 



where turnstiles are crucial for understanding large-scale diffusion [ 124 - 1271 ] . The escape region introduced above in 
the context of the plus-minus method represents precisely one half of such a turnstile. The main idea is to study 
the interaction of turnstiles, i.e., by varying the slope it shall be investigated whether one obtains "good" or "bad" 
conditions for particles to get from one turnstile into another, or maybe even to get mapped successively through a 
series of turnstiles. As before in case of the plus-minus method, such dynamical correlations are expected to show up 
in the curve for the parameter-dependent diffusion coefficient D(a). The advantage of turnstile dynamics is that it 
can be made quantitative by exemplifying all turnstiles with certain points of these regions. For instance, the peak 
of the turnstile one starts with may be represented by the critical point, x = 1/2. Now, one can try to compute the 
slopes for which this point maps into other turnstiles again, being exemplified by certain points, after certain numbers 
of iterations. 

This has been done in detail for the region 2 < a < 3, as shown in Fig. [l^. The dashed line in the figure represents 
the prcdiciton of D(a) for a simple random-walk model as suggested in Rcf. [46[, which is discussed in detail in Refs. 
p8| , [77|j78| . Note that, on a large scale, the model correctly accounts for the behavior of D{a) near a = 2, but that with 
respect to any fine structure, such a simple model is clearly totally apart. One can see three distinct series of values of 
a in the figure. To understand the nature of these series, one should consider the orbit of the critical point. The first 
iterate of x = 1/2 is in the second interval, (1,2). The series a values of a are defined by the condition that the second 
iterate of x = 1/2 is at the leftmost point of the upward turnstile in the second interval (1, 2) (a = 2.732), or that the 
third iterate is at the corresponding point in the third interval (a — 2.920), etc. The numbers on the D(a)-curve refer 
to the number of intervals the image of x — 1/2 has travelled before it gets to the appropriate point on the turnstiles. 
Series (3 points are defined in a similar way, but they are allowed to have two or more internal reflections within an 
interval before reaching the left edge of a turnstile. Series 7 points are defined by the condition that some image of 
x = 1/2 has reached the rightmost edge of an upward turnstile, i.e., some point x — k + 1/2, where k is an integer. 
One observes that each series produces a cascade of apparently self-similar regions of decreasing size, as the limits 
a — > 2 or a — > 3 are approached. These cascades provide a basis for a physical understanding of the features of D(a) 
in this region: Particles leave a particular unit interval through a turnstile and undergo a number of iterations before 
they are within another turnstile. Whether they continue to move in the same or in the reverse direction at the next 
and later turnstiles is a sensitive function of the slope of the map. Thus, the fractal structure of the D(a) curve is 
due to the effects of long-range correlations among turnstiles, and these correlations lead to changes of D(a) on an 
infinitely fine scale. We note that another way to understand this fractal structure is in terms of so-called "pruning" 
of the microscopic deterministic dynamics. That is, by varying the parameter a certain types of orbits may suddenly 
disappear. This means that with respect to a given symbolic dynamics of the map certain symbol sequences, which 
identify the orbits, do not exist anymore. This can be related to the irregularities of the diffusion coefficient |],|5(|. 

One should note that series 7 points completely label the maxima of higher order introduced before, and series (3 
points mark the respective minima. This way, in the region of the slope a < 3 the picture of quantitative turnstile 
dynamics is in full agreement with the results obtained by the qualitative plus-minus method outlined above (the 
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agreement has been checked to persist at least up to a level of extrema of second order). However, the application 
of turnstile dynamics has its limits: First, this method is of no use anymore for any higher value of the slope above 
a = 3. Thus, there is no other understanding of the structure in this range than the one provided qualitatively by the 
plus-minus approach. And second, even for values below a = 3 turnstile dynamics is quantitatively not completely 
correct: Apart from the lack of explaining the existence of the overhang above a — 3, a detailed analysis reveals 
further "tiny overhangs" at maxima of higher order, as, e.g., right above the maximum of first order in the region 



2<a<3ata = 2.414 [128 1 . In other words, the "turnstile values" marked in Fig. |lO|by series 7-points represent not 
the exact local maxima of higher order of the curve. The true local maxima are in fact shifted slightly to the right 
from these points, as in case of a = 3. The phenomenon of overhangs are further elucidated in Ref. J5^ , |l20[ . However, 
apart from the qualitative remarks in the context of the plus-minus approach and the additional insight provided by 
the approach in Ref. [5^J12C[ | , a detailed explanation of these overhang effects is still missing. 

Since turnstile dynamics points seem to separate self-similar regions, it is suggestive to use them as a tool to do some 
scaling. Series a-values are especially suitable for this purpose, because they form a scries of points which converges 
monotonically to a = 3, defining self-similar regions of decreasing size. These regions have been scaled according to 
the size of the intervals of the slope and of the respective diffusion coefficient intervals defined by 

Act; := a(series a-point (i + 1)) — a(series a-point i) 

ADi := D (series a-point (i + 1)) — D (series a-point i) , i = 0, . . . , 4 (40) 
(note that a — 2 has been taken here as series a point 0). To obtain scaling factors, the fractions 

SaW := — r and s D (i) := (41) 

Aa, AZ/; 

have been computed. They led to the two series of values given in Table Q. s a (i) seems to converge quite rapidly to a 
value around 3, whereas so{i) approaches not that fast a value apparently between 2 and 3, but this is no more than 
a guess based on the first five values of two infinite series, since the data set of slopes and D(a)-values is not sufficient 
for obtaining better results. At least these values suggest that some quantitative scaling is possible in this region. 

At this point, it should be stressed that the region below a = 3 is special, compared to any other region of the 
slope: Firstly, the structure of the curve is remarkably simple, as shown in Fig. [l0| Secondly, the number of wiggles 
of higher order is not constant with increasing order, but grows according to the structure described by the turnstile 
dynamics performed above. This is in contrast to the behavior of D(a) in the ranges 4 < a < 5 and 6 < a < 7, where 
one may have expected similar generalities. Thirdly, the region below a = 3 is the only one which is simple enough 
such that turnstile dynamics can successfully be applied at all, and this region seems to provide some simple scaling 
laws. All this nice behavior suddenly breaks down at the value a = 3, which is marked by the largest overhang of 
the whole curve. Therefore, it might be assumed that the point at a — 3 separates regions of fundamental different 
dynamical behavior of the map, i.e., the dynamics seems to be sufficiently simple below, but suddenly gets quite 
complicated above this value. In fact, there is further evidence that such a transition exists, as is discussed in detail 
in Refs. §§|r§]. 

Finally, another interesting feature of the parameter-dependent diffusion coefficient for map £ shall be pointed out. 
To obtain Fig. [ll], the first "derivative" of the D(a)-curve of Fig. [s] has been computed with respect to the full data 
set of D(a)-values available. This has been done in linear approximation, i.e., two adjacent points of the data set have 
been connected with straight lines, and the "derivative" D'(a) := AD / Aa for Aa <C 1 has been computed. However, 
since the curve of D(a) is expected to be nowhere differentiable with respect to the slope, the derivative defined 
above should be better denoted as a "pseudo-derivative", i.e., as a kind of mean value of AD over tiny regions of 
Aa which are approaching zero. Thus, the pseudo-derivative presented in Fig. [ll] is mathematically not well-defined. 
Nevertheless, the numerical results turn out to be quite reasonable in the sense that they reflect to a certain degree 
the fractal structure of the actual D(a)-curve. For instance, the three "bands" right above a = 2 seem to be due 
to the "triangle-like" self-similar structure of the region presented in Fig. [ll] and the "bursts" at odd integer values 
correspond to the occurrence of local maxima and to the deformations of simple self-similar structures. We remark 
that in Fig. [ll] the number of slopes for which D(a)-values have been calculated is not homogeneously distributed 
over the whole range 2 < a < 8, but that the number of points per interval is greater in certain regions of the slope, 
as, e.g., right above a = 2, a = 3 and around a — 5.6. These regions thus show up as slightly pronounced parts in the 
derivative plot, howe ver, t he structure of the curve does not seem to be affected by it. More details of this curve are 
discussed in Refs. [p||l20|. 



V. SUMMARY 
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A. Conclusions 



(1) A simple model for deterministic diffusion has been discussed where the microscopic scattering rules can be 
changed by varying a single control parameter. The diffusion coefficient of this model has been computed for a broad 
range of parameter values and shows a fractal structure as a function of the slope of the map. This result appears to 
be the first example of a dynamical system whose diffusion coefficient has an unambiguously fractal structure. 

(2) A general method to compute parameter-dependent diffusion coefficients for a whole class of piecewise linear 
maps has been developed, ft is based on the first passage method, which provides the definition of the diffusion 
coefficient for the dynamical system, in combination with the use of Markov partitions and transition matrices, 
which have been employed to solve the Frobenius-Perron equation of the dynamical system. For periodic boundary 
conditions, the parameter-dependent diffusion coefficient could be related to the second largest eigenvalue of the 
topological transition matrix. This method provides analytical solutions in simple cases and is also accessible to 
numerical implementations. 

(3) The method described above has also been applied to absorbing boundary conditions. Long-range boundary 
layers have been found in the eigenmodes of the deterministic dynamical system. They also show up in quantitative 
calculations of the diffusion coefficient. 

(4) Certain limits of the first passage method in combination with the use of transition matrices have been discussed: 
Drawbacks are especially the restriction to certain initial probability densities suitable for the application of transition 
matrices, as well as the "external" definition of the diffusion coefficient by the "matching eigenmodes" procedure of 
first passage. It turned out that this procedure is not well-defined anymore for smaller eigenmodes of the dynamical 
system. 

(5) A systematical way to find Markov partitions for the class of maps under consideration has been developed. 
This method has been used as the basis for computing the parameter-dependent diffusion coefficient for the dynamical 
system mentioned above. For this map, as well as for the whole class of maps under consideration, Markov partitions 
are conjectured to be dense in the set of parameter values. 

(6) A large number of eigenvalue problems of topological transition matrices, based on Markov partitions, has 
been solved numerically to compute the parameter-dependent diffusion coefficient for the model. In the course of 
these calculations, the reliability of well-known standard software routines for computing eigenvalue spectra has been 
checked critically, and numerical uncertainties have been pointed out. 

(7) Certain large- and small-scale structures in the eigenmodes of the topological transition matrices have been 
found. The large-scale structures support the existence of statistical diffusion in the dynamical system, whereas the 
small-scale structures refer to the specific microscopic deterministic dynamics of the model system. These results 
suggest that the strength of the fractal diffusion coefficient is related to the fine-scale structure of the eigenmodes. 

(8) A conjecture about the existence of diffusion coefficients for a broad class of one-dimensional maps has been 
made. This conjecture may shed more light on the origin of diffusion generated by a simple deterministic dynamical 
system and may show a way how to put the theory outlined in this article onto more solid mathematical grounds. 

(9) Qualitative explanations for the structure of the parameter-dependent diffusion coefficient over the full range 
of parameter values have, to a certain extent, been provided by simple heuristic considerations. 

(10) A more refined "turnstile dynamics" has been developed as a more quantitative approach to explain the 
structure of the parameter-dependent diffusion coefficient. It works in certain regions of the parameter values and 
provides a starting point for a scaling of certain self-similar structures. 

(11) By employing these qualitative and quantitative methods, certain interesting features of the diffusion coefficient 
have been discussed, i.e., the phenomenon of "overhangs" at local extrema, and the special simple character of an 
"initial region" for small parameter values, where diffusion sets in. 

(12) The numerically computed pseudo-derivative of the parameter-dependent diffusion coefficient seems to provide 
another characteristic property of fractal diffusion coefficients. 

B. Outlook 

The class of one-dimensional piecewise linear maps we have studied here by analyzing an example appears to be 
the most simple type of deterministic diffusive systems one can think of. Nevertheless, we have shown that the 
diffusion coefficient of such a map changes in a fractal way by varying a control parameter. Starting from this 
fundamental result, there are at least two directions in which our research can be pursued: One way is to study 
whether other transport quantities, like electric conductivities, chemical reaction rates, or magnetoresistancies can 
exhibit such an irregular behavior as well. Another way is to investigate whether fractal transport coefficients exist 
in more complicated, and thus more realistic, dynamical systems. 



16 



First steps in these directions have already been taken: For example, a bias has been added to the simple map 
discussed in this paper. This g enerates an average current of particles which again exhibits fractal structures by 
varying t he bias as a parameter [104,121]. Moreover, for small enough bias the current can run oppo site to the bias 
[104, 121 1, and for other parameter values the diffusion coefficient is zero with nonzero current [ 104 1 . Deterministic 
diffusion coefficients and deterministic currents which change irregularly by varying respective parameters have also 
been found in parameter-dependent two-dimensional multibaker models In their transport properties, these 

models are closely related to the class of one-dimen siona l maps discussed here. However, in addition they provide 
more physical features like being time-reversible |3£,104] and being area-preserving or dissipative in a well-defined 
sense. Negative currents have been observed in these systems as well, and a parameter-dependent diffusive-reactive 
multibaker yields chemical reaction rates which are also fractal in a parameter |3Sj| . 

In the periodic Lorentz gas, however, there is up to now no clear indication about a fractal behavior of transport 
coefficients. First results of computer simulations for the diffusion coefficient with respect to varying the density of 
scatterers show a very smooth c urve , which indicates that if there are fractal fluctuations in the parameter at all they 
must occur on a very fine scale 12£]. On the other hand, for the thermostatted periodic Lorentz gas with an electric 



field computer simulations of several groups show clearly a very irregular behavior of the conductivity by varying the 
field strength|2(],H,H,|l}. 

In one case, the numerical results could even be confirmed by calculations based on cycle 

expansions ||31|. 

Whether deterministic phenomena of this kind play a role in real statistical, ex periment ally accessible systems 
is a very open question. Following the chaotic hypothesis of Gallavotti and Cohen [ 130 ,131 1 , one may believe that 
these phenomena are rather due to the simplicity of the models and should eventually disappear if the systems are 
getting more complex. Respectively, we would expect that certain necessary conditions must be fulfilled for systems to 
exhibit characteristics of fractal transport coefficients which may contradict the spirit of the chaotic hypothesis, such 
as being spatially periodic, low-dimensional, and such that particle-particle interactions are not of main importance. 
Physical systems of this kind could - to a certain extent - already be realized experimentally in form of so-called 
antidot lattices. Here, magnetoresistancies which fluc tuate irregularly by varying the field strength have already been 
observed experimentally in a classical limit [132,133], and t o a certa in extent they have been explained theoretically 
by identifying special orbits in the microscopic dynamics [134,135]. Another candidate of a system where certain 
irregularities in transport could be of a deterministic origin are so-called ratchets, where negative currents have 
already been found experimentally, as well as in theoretical models (see, e.g., Refs. [136,137] and further references 
therein). In fact, it can be argued that there exists a relation between certain types of ratchets and the class of 
one-dimensional maps (supplemented by a bias) studied here [ 104, 138 1. 

Fractal transport coefficients in one-dimensional maps actually appear to be stable with respect to imposing dif- 
ferent kinds of random perturbations on the system [104]. This means that the fractal structure gradually smoothes 
out by increasing the perturbation strength and thus survives in form of irregular oscillations on finer scales if the 
perturbation is small enough. Although there are exceptions to this behavior [139], we believe this to be the typical 
scenario of how a possible fractality of parameter-dependent transport coefficients may appear if the system is not 
completely deterministic. 
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APPENDIX A: TRANSITION MATRIX METHOD FOR CALCULATING DIFFUSION COEFFICIENTS 



In this Appendix, we deal with the problem of how to solve analytically the eigenvalue problems of the topological 
transition matrices T(a) formally introduced in Section III A. These transition matrices are the key ingredients of 



the Frobenius-Perron matrix equation, Eq. (^), and their eigenmodes and eigenvalues determine the deterministic 
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diffusive dynamics of the map at the respective parameter value of the slope, as outlined in Section |l| Two simple 
examples of such transition matrices have already been given by Eq. ([TJ) an d Eq. (|24|), the problem of transition 
matrices for more complicated Markov partitions has been discussed in Section III B . 

In the first subsection of this Appendix, we will consider integer slopes with periodic boundary conditions. In this 
case, trivial Markov partitions can be found for which the corresponding transition matrices have a simple structure. 
We analytically solve the eigenvalue problems of these transition matrices by standard methods and compute the 
parameter-dependent diffusion coefficients as defined by the first passage method. In the second subsection, we solve 
the eigenvalue problems of the respective transition matrices for two integer slopes with absorbing boundary conditions, 
namely for a = 3 and a = 4. In the third subsection, we discuss the problem of non-trivial Markov partitions for non- 
integer slopes. We depict Markov partitions for some special values of the slope, sketch the corresponding transition 
matrices, outline how to solve their eigenvalue problems analytically, and give the results for their eigenvalues and 
diffusion coefficients. 



1. Integer slopes with periodic boundary conditions 



As discussed in Section [II, for periodic boundary conditions all topological transition matrices corresponding to 
Markov partitions are block circulants. For certain values of the slope it is possible to r e duce these block circulants 
to simple circulants for which it is known how to solve their eigenvalue problems [100 102 1. We employ here the 
approach of Berlin and Kac as described in App. A of Ref. 5 [10C]. We summarize their main formulas and use them 
first to solve the problem of even integer slopes for periodic boundary conditions. We then treat analogously the 
slightly more complicated case of odd integer slopes. 

Let T be a cyclic matrix of the type 



T := 
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We want to calculate the eigenvalues Xm and the eigenvectors \4> m > associated with T, that is, 



T \<j) m >= Xrn \4>m > 
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(A2) 



According to Ref. [10C], the eigenvalues are given by 
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and the corresponding eigenvectors are 
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(A5) 



with a m and b m to be fixed by suitable normalization conditions. 

We now compute the diffusion coefficient for all even integer slopes a — 2k , k € N, of map L by solving the 
eigenvalue problem of the respecti ve gen eral transition matrix T(a) of the system. This transition matrix can be 
constructed as discussed in Section III A by using the same Markov partition for all slopes, that is, the one which is 
depicted in Fig. We find that the matrix corresponding to this Markov partition is a simple circulant with matrix 
elements of 
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h = 2 , t 2 = 1 , 



ta/2 — 1 ! t a /2+l — . 



! ii-a/2+1 — , ^L-a/2+2 — 1 j • • • i = 1 



(A6) 



The (s + a/2 — l)th row, 1 < s < L — a/2, of the corresponding eigenvalue problem defined by Eq. (A2) is thus 
determined by 



+ . . . + 4>T ,2 ~ 2 + ^T' 2 - 1 + 4ti a/2 + ■■■ + 4>T~ 2 = &C /2_1 



According to Eq. (A3), the eigenvalues are 

Xm = 2 + r m +... + r<T 2 )/ 2 + r^(- 2 )/ 2 + . . . + r 



L-l 
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and by using Eq. (A4) we obtain explicitly 
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The corresponding eigenvectors are given by Eq. (A5). Note that the largest eigenvalue xo is equal to the slope of the 
map. This is a consequence of the fact that the topological transition matrices discussed here can be map ped o nto 
stochastic transition matrices and that the Perron- Frobenius theorem for non- negative matrices applies |32,140|. It 
can be proven that this property holds for any topologic al tra nsition matrix of map C which is defined on the basis 
of Markov partitions with periodic boundary conditions [104]. According to the matrix Frobenius-Perron equation, 



Eq. (|l^), the corresponding largest eigenmode determines the equilibrium state of the system, which is here simply 
uniform. In analogy to Eq. (|2l]), the second largest eigenvalue gives the decay rate, and the respective second largest 
eigenmode governs the diffusive transport in the map. For slope a = 2k we thus obtain a decay rate of 



ldec(a) = In 
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~T? 24 



With Eq. (9), this gives a diffusion coeficient of 
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We now want to do the same calculation for all odd integer slopes a — 2k — 1 , k E N, again by employing periodic 
boundary conditions. As a Markov partition for all these slopes we use the same partitioning underlying the transition 
matrix as given by Eq. (p4]), that is, its parts are all of length 1/2. Here, the corresponding a-dependent transition 
matrix is a block circulant where the single blocks consist of 2x2-matrices. To illustrate the general structure of this 
matrix, we give the first rows and columns of the special case a — 5 which is 



T(5) 



/1101000000 
1 1 1 1 
1 1 1 1 
10 110 10 10 
1010110100 
0010110101 
0010101101 

V i i : 
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We now write down the eigenvalue equation, Eq. ( |A.2j ), for this matrix. By using the notation itj^ for the fcth odd 
component of \<fi m > and for the feth even component we obtain for k odd, 1 < k < L — 3 , < m < L — 1, 



fc+1 , k+2 , fc+2 , fc+3 _ fc+2 
'm ^ u m T V m T u rn — XmU m 



(A13) 



supplemented by periodic boundary conditions for the respective first and last rows of the matrix. This leads to 



,fc — ^.fc+i 



, yielding an equation for wL which reads 
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(A14) 
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again by providing respective peri odic boundary conditions. Thus, we have reduced the eigenvalue prob lem for the 
initial block circulant of Eq. ( |A1S| ) to the eigenvalue problem of a simple circulant as given by Eq. ( A14 ). The same 
reduction procedure can be carried out for general odd integer value of the slope. The reduced eigenvalue equation 
then reads 

U f n + ...+^„ +a - 1 = Xm ^ +(a - 1)/2 , (A15) 
supplemented by respective periodic boundary conditions. With the Berlin-Kac method we obtain for the eigenvalues 

(o-l)/2 



Xm = l + 2 J2 cos(fl m s)~a — 1 24 ' (L^oo) , 0<m<L-l . (A16) 
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The eigenvectors are given here by 
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L 

with a rn and b rn to be fixed by suitable normalization conditions. For the decay rate we get 

7dec (a)c^^^ (L^oo) , (A18) 



which leads to a diffusion coeficient of 



£(<*) = . (A19) 



2. Integer slopes with absorbing boundary conditions 

For absorbing boundary conditions, the corresponding transition matrices are no block circulants anymore, but 



they belong to the broader class of banded square block Toeplitz matrices, as pointed out in Section III A. For these 
matrices no general methods are known for solving their eigenvalue problems analytically. However, in certain cases 
analytical solutions can still be obtained by straightforward calculations, as we will show for the two integer values of 
slope a = 3 and a = 4. 

We first consider the case a — A. The transition matrix for this parameter value is identical to the one given by 
Eq. (S3) except that the upper right and the lower left corners are filled with zeroes because of absorbing boundaries. 
The eigenvalue problem of this matrix can now be solved in analogy to the calculations performed by Gaspard in Ref. 



52[. The eigenvalue equation, Eq. (A2), reads here 



4> L m + 2rtt 1 +€: 2 =Xrn€t 1 . 0<fc<L-l , (A20) 

supplemented by the absorbing boundary conditions <^> m = 0^ +1 = 0. Since this equation is of the form of a discretized 
ordinary differential equation of degree two we make the ansatz 

</4 = a cos(fc6») + b sin(fc0) , < k < L + 1 . (A21) 

The two boundary conditions then lead to 

a = and sin((L + 1)9) = (A22) 



yielding 



L + l 

The eigenvectors are then determined by 
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m = T—r , l<m<L . (A23) 
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(A24) 



with b as the normalization constant. Putting this equation into Eq. ( |A20| ) gives as the eigenvalues 

Xrn = 2 + 2 COS d m . 



(A25) 



Note that in case of absorbing boundary conditions the largest eigenvalue is not equal to the slope of the map, but 
determines the escape rate of the system. Correspondingly, the largest eigenmode is the diffusive mode of the map. 
However, for L — > oo the largest eigenvalue goes to the exact value of the slope, which therefore serves as an upper 
limit of the eigenvalue spectrum. This is conjectured to be true for any topological transitio n m atrix of map C 
which is defined on the basis of Markov partitions by employing absorbing boundary conditions [104]. In th e limit of 
chainlength L — > oo these results lead to the escape rate and diffusion coefficient presented in Section III A, Eqs. ( |30| 

We now treat analogously the case of slope a = 3 for absorbing boundary conditions. We know from the previous 
subsection that for odd integer slopes the Markov partitions, and thus the respective transition matrices, are a bit 
more complicated. The transition matrix for a = 3 is identical to the one given by Eq. (^4|) except that the upper 
right and lower left corners are filled with zeroes because of absorbing boundaries, as before. 

To write down the eigenvalue equation, Eq. (A2), for this matrix we use the same notation as in the previous 



subsection for odd integer slopes. With u 1 ^ being the fcth odd component of \<j> m > and v m being the fcth even 
component we obtain for k odd, < k < L, 



ul + u 



k+l 



fe+1 
m 
,,fc+l 



XmU 



k+l 



supplemented by the absorbing boundary conditions 
equation for wL which reads 



,L+1 



(A26) 

0. This again leads to = u^ 1 , yielding an 



,fc+i 



,k+2 



XmU 



k + l 



< k < L 



(A27) 



with the respective absorbing boundary conditions. We again use Eq. (A.21) as an ansatz to solve this equation. The 
two boundary conditions then lead to 



and sin((Z + 2)0) = 



(A28) 



yielding 



1 < m < L + 1 



(A29) 



The eigenvectors are then determined by 



u m = bsm(k9 m ) 

<4 = &sin((fc + l)0 m ) 



(A30) 



with b as the normalization constant. Putting this equation into Eq. (A27) gives as eigenvalues 

X m = l + 2cos0 m . (A31) 
In the limit of chainlength L — > oo this leads to the results for escape rate and diffusion coefficient presented in Section 



III A, Eqs. (30 -B2J) 



We remark that we do not have analytical solutions of the eigenval ue pro blems for integer values of the slope above 
a = 4 with absorbing boundary conditions. Here, the ansatz of Eq. ( |A2l| ) does not seem to be sufficient because of 
long- range boundary layers which are induced by the absorbing boundary conditions, see also the remarks in Section 
III A to this problem. 



3. Non-trivial Markov partitions with periodic boundary conditions 

In this subsection, we discuss some examples of non-trivial Markov partitions where analytical solutions of the 
respective eigenvalue problems of the transition matrices can still be obtained analytically. This can be done by 
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reducing block circulants onto simple circulants in the way illustrated in the previous subsections for odd integer 
slopes. In the following, we will only consider periodic boundary conditions, because then the general method of 
Berlin and Kac to solve the respective eigenvalue problems can be applied. 

We discuss two different scries of Markov partitions. For the simplest case of the first series we briefly outline of 
how to perform the calculations, and we give the results for the eigenvalues and the diffusion coefficient. For the next 
parameter value of this series we give only the main results, before writing down the respective general formulas for 
the whole series. For the second series we only deal with the first two parameter values by giving the main results. 

In each case we proceed by first depicting a box map of the full chain of boxes with its Markov partition in a figure. 
We indicate of how the respective Markov partition has been computed and give the exact value of the slope by which 
it is defined. On this basis, we sketch the corresponding transition matrix for the full chain of boxes and give the 
main results for eigenvalues and diffusion coefficient. 



a. Series 1 



As has been pointed out in Section III B , a Markov partition is defined via a generating orbit which obeys Eq. 
(|36|). By using the notation of this equation, the first series of Markov partitions discussed here is characterized by 
6 = and the number of iterations of the reduced map being n = 1. The first case of this series is obtained from the 
solution of this equation for the slope a being restricted between 2 and 4. The second case refers to 4 < a < 6, the 
general case is for solutions 2k < a < 2(k + 1) , k € N. 



Case 1 : 

This is the simplest case and corresponds to the smallest slope of this series, as illustrated in Fig. |l2| (a). As one 
can infer from the figure, the precise value of the slope can be computed from the equation 



l = 2(l + e)e 0<e<l/2 
where a = 2(1 + e), which leads to the solution 

a= (V3-l)/2~ 2.73205 



(A32) 



(A33) 



The partition of the full chain of boxes can be constructed by continuing the box map of Fig. [1^ ( a) per iodically. The 
corresponding transition matrix can then be obtained from this partition as described in Section [II B and reads 



T(2. 73205) 



(A34) 



/l 1 

1 1 

1 1 1 

1 1 1 
0001010000 

1 1 1 

1 1 1 

V; ; i j 

where the first, second, third, . . ., three rows correspond to the first, second, third, . . ., box of the chain. We now 
reduce this block circulant to a simple circulant. Since the Markov partition of each box consists of three parts, we 
use t hre e different symbols uj^ , , as components of the eigenvectors to write down the eigenvalue equation, 
Eq. (|a|), of this matrix, 



,k+l 



to. 



fc+1 
m 
fc+1 



Xm v rn 
fa 

XmW ri 



Xm^ 



fc+1 



1 < k < L 



(A35) 



supplemented by respective periodic boundary conditions for the first and the last row of the matrix. From these 
equations it is immediately obtained that 



^m +1 + X m V k m 



,fc+l _ 



"m + XmV^ 1 



(A36) 
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which leads to 



2v k +v k ~ 1 = y 2 v k 



m 1 m 



(A37) 



with respective periodic boundary conditions. This again defines the eigenvalue problem of a simple circulant which 
we can solve by the methods used before. For the eigenvalues we obtain 



Xm = 1 ± y/l + 2 cos0 m ~ 1 ± V3(l - : f) (L-oo) , 

6 

where the second largest eigenvalue yields a diffusion coefficient of 

V3 



2nm 



< m < L - 1 



L>(2.73205) 



6(1 + V3) 



0.10566 



(A38) 



(A39) 



As pointed out above, the largest eigenvalue is again identical to the slope of the map. It is related to an equilibrium 
eigenmode which is here a periodically continued piecewise constant function, based on the single parts of the Markov 
partition. 

Case 2 : 

The second case of this series is the Markov partition defined by the respective value of the slope between 4 and 6. 
Its box map partition is illustrated in Fig. 12 (b) and corresponds to a = 2 + y8 ~ 4.82843. The transition matrix 
reads 

\ 



T(4.8284) = 
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(A40) 



/ 

Again, this block circulant can be reduced to a simple circulant. By some calculations which are quite analogous to 
the ones of the previous example we obtain as eigenvalues 



Xm = 1 + cos m ± a/(1 + cos 8 m ) 2 + 2 cos(2# m ) + 2 cos(3# m ) 
9 m and m as before, with a diffusion coefficient of 



2 - 
2 - JUL ± ^8(1 
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oo 



D(4.8284) = 



fV2 



4(l + x/2) 



0.65273 



(A41) 



(A42) 



General case of this series : 

We now give the general solution for eigenvalues and diffusion coefficient of all cases of this series of Markov 
partitions. Let 



1 



a(p) = 2(p + 1 + e) with e = - 1 + v / 2p + 2 + (p+ l) 2 ) , p e N 



(A43) 



be the slope of map C. Then a(0) is the value of the slope of case 1 and a(l) is the value of the slope of case 2. Any 
higher value of a(p) , p > 1 , defines a Markov partition which is of the same type as in the two previous examples, 
that is, it fulfills the general definition of this series given at the beginning. For this series of Markov partitions the 
general eigenvalue problem of the corresponding transition matrices can be solved by generalizing the calculations 
above. This yields for the eigenvalues 
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Xm = 1 + ^ cos(fc6» m ) ± 



fc=l 



\ 



p+1 



[(l + ^cos(fc6» m )) 2 + 2^cos((p + /c)6» m ) , 0<m<L-l 



(A44) 



fc=i fc=i 
In the limit of L — > oo, this leads to a diffusion coefficient of 

{2p 2 + p) v / (l+p){3+p) + 2p 3 + 17p 2 + 20p + 6 



£>(p) 



12(^(1 +p)(3+p) + 3 +p) 



(A45) 



6. Series 2 



By again referring to Eq. (p6|), and using the notation of Section |III BJ , a second series of Markov partitions is 
defined by 8 = 1 — e and the number of iterations of the reduced map being n = 1 . In the following we give the main 
results for only the first two cases of this series, which are based on the solution of this equation for the slope a being 
2 < a < 4, and 4 < a < 6, respectively. 



Case 1 : 

The partition illustrated in Fig. O (c) corresponds to the slope 



a = 1 a 2.56155 



(A46) 



and defines, periodically continued, a transition matrix of 



T(2.56155) = 
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(A47) 



The eigenvalues are 



1 /9 „ 

Xm = - ± A / ^ + 2 COS I 



2-7rm 

'm , 9 m = — , < m < L - I 



(A48) 



the diffusion coefficient is 



L>(2.56155) 



171+V17 



0.09468 . 



(A49) 



Case 2 : 

The partition illustrated in Fig. |l2| (d) corresponds to the slope 



a =1±£U 4.70156 



(A50) 



and defines, periodically continued, a transition matrix of 
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T(4.70156) 
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The eigenvalues are 



l + 2cos6> m ± /l^ + 2cos6/m ^ 2 + 2 + 2cos(9 m + 2cos(26» m ) + 2cos(30 m ) , 



9 m and m as before, the diffusion coefficient is 



£(4.70516) = 



f + 



31 

vTT 



0.62f22 



(A52) 



(A53) 



We finally remark that analytical calculations which are similar to the ones performed here have been carried out in 
Refs. |6l pC| , |5f|| . These calculations are based on cycle expansions, and the diffusion coefficient has been computed for 
different piecewise linear maps. 
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1 


2 


3 


4 


Sa(i) 


3.902 


3.423 


3.186 


3.079 


3.033 


sd(i) 


1.128 


1.510 


1.721 


1.892 


2.038 



TABLE I. Scaling factors for the initial region 2 < a < 3 of the parameter-dependent diffusion coefficient D(a) (see text). 
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FIG. 1. Sketch of a simple one-dimensional random walk model on the real line, in contrast to deterministic dynamics. 




FIG. 2. Illustration of a simple model of deterministic diffusion, see the dynamical system map C, Eqs.(g) to (j6|), for the 
particular slope a — 3. The dashed line refers to the orbit of a moving particle. Its initial condition is indicated by a black 
arrow close to the x-axis. The particle moves under the action of the one-dimensional piecewise linear map shown in the figure 
by jumping from box to box. 
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FIG. 3. Partition of map C at slope a — 4 (dashed grid). 
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FIG. 4. The two second largest eigenmodes of map chainlength L — 100, for slope a — 3 with periodic boundary conditions. 
The eigenmodes exhibit a step-like fine structure and differ by a phase shift. In comparison, the respective two eigenmodes 
obtained from solving the diffusion equation Eq.(Q) have been included as dashed lines. They are almost indistinguishable from 
the map eigenmodes. 
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FIG. 5. Largest eigenmodes of map C for odd integer values of the slope a with absorbing boundary conditions and comparison 
to the largest eigenmode of the diffusion equation Eq.(W). The inset is a magnification of the region around x = 0. 
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FIG. 6. Four examples of non-trivial Markov partitions of map C at different values of the slope. Diagram (a) is for the 
slope a ~ 2.057, (b) for a ~ 2.648, (c) for a ~ 6.158 and (d) for a — 7.641. In (a) and (b) the bold black lines with the arrows 
show the generating orbits of the partitions, that is, the orbits which define the single partition points (see text). The two large 
arrows below the orbits indicate the respective initial positions for the generating orbits. 
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FIG. 7. Second largest eigenmodes of map £ at two non-trivial Markov partition values of the slope with periodic boundary 
conditions: full modes and magnifications of their fine structures. For both parameter values there are two largest eigenmodes 
which differ by a phase shift. Diagrams (a) and (b) are for slope a — 3.0027, chainlength L = 100, (c) and (d) for a = 2.0148, 
chainlength L = 90. 
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FIG. 9. First qualitative approach to understand the structure of the parameter-dependent diffusion coefficient D(a), denoted 
as plus-minus method (see text). The variation of the microscopic scattering process via changing the slope a by Aa is 
heuristically related to variations in the strength of the diffusion coefficient. The plus (minus) signs refer to subintervals 
where forward (backward) scattering occurs at the next iteration of the map (particle moves to the right, or left, respectively). 
The qualitative argument is that the sequence of dominant forward or backward scattering by varying the parameter induces 
oscillations in the strength of the diffusion coefficient. Figs, (a) and (b) are for one iteration, Figs, (c) and (d) for two. 
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FIG. 10. Blowup of the region of slope a < 3 for map C with the solution for a simple random walk model (dashed 
line) and labels for parameter values which are significant for "turnstile dynamics" (see text). Turnstile dynamics establishes a 
quantitative relation between the local maxima and minima of the parameter-dependent diffusion coefficient and the underlying 
microscopic chaotic scattering process. For some parameter values, the turnstile coupling is shown by pairs of boxes. The graph 
consists of 979 single data points. 
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FIG. 11. Numerical derivative of the parameter-dependent diffusion coefficient D(a) of map C with respect to the slope a. 
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